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Alien species are of major concern to conservation biology, 
agriculture and the human society, as they may become in¬ 
vasive and successfully compete with native species, nega¬ 
tively influence crop production, and cause health prob¬ 
lems. Several reptile species are considered invasive alien 
species, but most of them are little studied (Kraus 2009). 
Amongst these species are a few chameleons. Apparently, 
several of the Chamaeleo chamaeleo (Linnaeus, 1758) pop¬ 
ulations in the Mediterranean (e.g. Crete, Malta, Sicily) are 
non-autochthonous and their presence can be linked to 
historical trade routes (e.g. Klaver 1981). Chamaeleo af- 
ricanus Laurenti, 1768 from northern Africa was translo¬ 
cated from the Nile delta onto the Peloponesos where it has 
maintained reproducing populations perhaps since antiq¬ 
uity (e.g. Bohme et al. 1998, Bohme 2000). More recently, 
the Malagasy Furcifer pardalis (Cuvier, 1829) has become 
invasive on Reunion Island (e.g. Bourgat 1970). 

One of the most successful invasive chameleons is Jack- 
sons chameleon, Trioceros jacksonii (Boulenger, 1896), 
which is native to Afromontane Kenya and Tanzania 
(Fig. 1). Popular as a pet and commonly traded and im¬ 
ported by the thousands from Kenya and Tanzania (Car¬ 
penter et al. 2004), this species was introduced to Oahu, 
Hawaiian Islands in 1972 (e.g. McKeown 1991, Waring 
1997). Apparently, established populations are referable 
to the large, up to 35 cm in total length, subspecies T. j. 
xantholophus (Eason, Ferguson & Hebrard, 1988) from 
Kenya and Tanzania (McKeown 1991) (Fig. 2). As a result 
of inter-island transport which has only been restricted 
since 1997, Jacksons chameleon currently has established 
wild populations on all the main islands including the is¬ 
land of Hawaii, Maui and Kauii as a result of multiple intro¬ 
ductions and subsequent spread (Fig. 1; Waring 1997 cited 
in Hagey et al. 2010). The species is especially common in 
wetter habitats at elevations from 100-1,000 m above sea 
level (McKeown 1996, Holland et al. 2009) and inhab¬ 
its disturbed habitats as well as forested areas (McKeown 
1991, Waring 1997). The ecological impact of Hawaiian T. 
jacksonii has not been studied. Only recently, Holland et 


al. (2009) documented predation on endemic snail and in¬ 
sect species, including some threatened with extinction, 
and Hagey et al. (2010) commented on its foraging strat¬ 
egy. Further studies on the natural history and distribution 
of T. jacksonii on Hawaiian Islands, coupled with monitor¬ 
ing and management efforts, are required. 

GIS-based Species Distribution Models (SDMs) are a 
helpful tool towards these goals, as they focus on species 
potential distributions. In the case of invasive alien spe¬ 
cies, they can be understood as spatial risk assessments. 
Technically, a SDM assesses a species’ ecological niche and 
projects it into geographical space (Rodder et al. 2010a). 
Niche information is obtained at presence and sometimes 
absence sites of the study species; climatic niche informa¬ 
tion has been demonstrated to be a good predictor for in¬ 
vasion success in alien reptile species (e.g. Bomford et 
al. 2009; Rodder et al. 2009b). The resulting map shows 
similarities with the niche elsewhere, e.g. by different grid 
values, interpreted as ‘likelihood’ of the species’ potential 
occurrence (e.g. Peterson & Vieglais 2001, Rodder & 
Lotters 2010). Projecting SDMs onto past or future cli¬ 
mate scenarios may allow simulations of a species’ poten¬ 
tial response to environmental changes (e.g. Rodder 8c 
Schulte 2010). For example, Rodder (2009) provided a 
spatial risk assessment for the invasive non-autochthonous 
frog Eleutherodactylus coqui on the Hawaiian Islands, sug¬ 
gesting range allocation towards higher elevations under 
future anthropogenic climate change, thus affecting pro¬ 
tected areas. 

In this paper, we develop SDMs for Jackson’s chameleon 
in the Hawaiian Islands under current and future climates 
based on the species’ climatic envelope at its invaded and 
native African distribution ranges. The goal is to appreci¬ 
ate the geographical potential of this species in its invaded 
range. 

For SDM building, we used presence data from the na¬ 
tive African and the invaded Hawaiian ranges of Trioceros 
jacksonii as available via the Global Biodiversity Informa¬ 
tion Facility (http://www.gbif.org) and HerpNet (http:// 
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Figure 1. Species records used for model building in T. jacksonii s invaded Hawaiian range and its native range in Kenya and Tanzania 
(insert). 


www.herpnet.org). When necessary, georeferencing of lo¬ 
cality data was conducted with the BioGeoMancer (http:// 
bg.berkeley.edu/latest/). DIVA-GIS 7.1.6 (Hijmans et al. 
2002; http://www.diva-gis.org) was used to check accuracy 



Figure 2. Non-collected male Trioceros jacksonii xantholophus 
from Parklands in Nairobi. 


of coordinates. Records were only considered if they could 
be unambiguously assigned to a single grid cell, leaving a 
total number of 81 presence records for modelling (Fig. 1). 

Information on current climate, i.e., for the period 
1950-2000, was obtained from the WorldClim database 
at a spatial resolution of 30 arc sec (Hijmans et al. 2005; 
http://www.worldclim.org). For future climate, as expected 
for the year 2080, we downloaded interpolations of three 
different global change scenarios of the third and fourth 
assessments of the IPCC (2007); spatial resolution 30 arc 
sec (Ramirez & Jarvis 2008; http://gisweb.ciat.cgiar.org/ 
GCMPage). These anthropogenic climate change pre¬ 
dictions, simulating conditions as expected in the two 
IPCC story lines A2a and B2a, were derived from simu¬ 
lations of three climate models CCCma-CGCM2 (Flato 
et al. 2000, Flato & Boer 2001), CSIRO-MK2 (Gordon 
& O'Farrell 1997), and UKMO-HadCM3 (Gordon et 
al. 2000, Pope et al. 2000). Climate variables comprised 
monthly minimum and maximum temperatures, and aver¬ 
age precipitation per month. Based on these, we computed 
19 ‘bioclimatic’ variables with DIVA-GIS; these are com¬ 
monly used for species distribution modelling (e.g. Busby 
1991, Beaumont et al. 2005). Since multicolinearity of pre¬ 
dictor variables may hamper successful SDM projection 
through space or time (e.g. Heikkinen et al. 2006), we 
extracted all ‘bioclimatic’ values from the species records 
and performed a pair-wise correlation analysis based on 
Pearsons correlation coefficient with XLSTAT 2010 (http:// 
www.adinsoft.com). For SDM computation, only variables 
with R 2 < 0.75 were considered, resulting in a final variable 
set comprising ’mean monthly temperature range’, ’tem¬ 
perature seasonality’, ’maximum temperature of the warm¬ 
est month’, ’minimum temperature of the coldest month’, ’ 
annual temperature range’, ’annual precipitation, ’precipita¬ 
tion in the wettest month’, ’precipitation seasonality’, ’pre¬ 
cipitation in the warmest quarter’ and ’precipitation in the 
coldest quarter’. It was suggested that these variables de¬ 
scribe environmental conditions that are physiologically 
important to T. jacksonii (Bennett 2004, Andrews 2008, 
Lin 1979, Lin & Nelson 1980). 
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Figure 3. Comparisons of the climate spaces occupied in T. jack- 
soniis native and invaded Hawaiian ranges. 


For SDM development, Maxent 3.2.2 was applied (Phil¬ 
lips et al. 2006, Phillips & Dudik 2008), a machine¬ 
learning algorithm following the principle of maximum 
entropy (Jaynes 1957). Maxent has been successfully ap¬ 
plied in studies assessing possible spreads of invasive al¬ 
ien species (Giovanelli et al. 2007, Rodder & Lotters 
2010), pathogens (Rodder et al. 2009a, Lotters et al. 2010 


“2009”), and potential impact of climate change (Rodder 
& Weinsheimer 2009, Rodder et al. 2010b). In most cas¬ 
es, Maxent outperforms comparable SDM methods, even 
when the sample size is small (Hernandez et al. 2006, 
Wisz et al. 2008). This was further confirmed by Maxents 
ability to predict previously unknown populations in poor¬ 
ly known species (Pearson et al. 2007, Weinsheimer et 
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Figure 4. Potential distribution of Jacksons chameleon on the Hawaiian Islands under current climate and two future anthropogenic 
climate change scenarios (A2a, B2a), based on a Maxent Species Distribution Model and bioclimatic conditions at the species’ native 
African occurrence. Warmer colours indicated higher suitability to the species when considering climate conditions only (each left) 
and when considering climate as well as land cover (each right). Established reserves are indicated as crosshatched areas. 


al. 2010) and even species that are likely to be new to sci¬ 
ence (Raxworthy et al. 2003). Next to species presence 
records, Maxent requires a set of pseudo-absence records 
for SDM computation. The selection of these is a crucial 
step towards successful model building (Phillips 2008, 
Mateo et al. 2010, VanDerWal et al. 2009). Herein, we 
randomly selected 20 pseudo-absence records per species 
record within a circular buffer of 10 km as recommended 
by Mateo et al. (2010). This results in a set of pseudo-ab¬ 
sence records with a similar spatial structure as present in 


the species records, equalling the target group background 
approach suggested by Phillips et al. (2009). 

For model evaluation, 70 % of the species records were 
used for model training and the remaining 30 % for model 
testing via the Area Under Curve (AUC), referring to the 
Receiver Operating characteristic Curve (Phillips et al. 
2006, Fielding & Bell 1997). Furthermore, Maxent allows 
tracing back the relative contribution of each variable to the 
final model, and so facilitates comparisons with physiologi¬ 
cal properties of the target species. As data splits may affect 
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the model due to the exclusion of some species records that 
represent features not present in the remaining records, this 
procedure was repeated 100 times and averages were com¬ 
puted subsequently. Land cover data for the Hawaii Islands 
were obtained from the Coastal Service Center of the Na¬ 
tional Oceanic and Atmospheric Administration (http:// 
www.csc.noaa.gov/crs/lca/hawaii.html); land cover classes 
included high- and low-intensity development, cultivated 
land, grassland, evergreen forest, scrub/shrub, palustrine 
emergent wetlands, unconsolidated shoreline, bare land, 
and water. Information on current reserve networks was 
adopted from the Office of Planning, State of Hawaii. 

Comparisons of the occupied climate space in T. jackso¬ 
nii s native and invaded Hawaiian ranges revealed that the 
highest similarities existed in the precipitation in the wet¬ 
test month’ and the precipitation in the warmest quarter’ 
and the greatest differences in the mean monthly tempera¬ 
ture range’, annual temperature range’ and the precipita¬ 
tion in the coldest quarter’ (Fig. 3). The minimum tem¬ 
perature of the coldest month’ experienced within T. jack¬ 
sonii s invaded range is well above its critical thermal mini¬ 
mum (5.3 ± 0.5 °C, n = 8, Bennett 2004), but well below 
its voluntary minimum temperature (29.1 ± 0.32 °C, n = 
7, Bennett 2004). Its preferred body temperature (32.1 ± 
0.24 0 C, n = 7, Bennett 2004), voluntary maximum body 
temperature (34.2 ± 0.31, n = 7, Bennett 2004), panting 
threshold (36.2 ± 0.67, n =5, Bennett 2004), and critical 
thermal maximum (41.0 ± 0.15, n = 8, Bennett 2004; 42.0 
± 0.91 0 C, n = 10, Lin 1980) are slightly higher than the 
maximum temperature of the warmest month within its 
invaded range. Differences of up to 6 °C may be compen¬ 
sated by heliothermic thermoregulation, and feeding ac¬ 
tivity may occur at temperatures as low as 10 °C (Bennett 
2004), however. 

It has been suggested that climatic conditions support¬ 
ing successful reproduction and balanced sex ratios can be 
strong predictors for the geographic range limits of rep¬ 
tile species (e.g. Rodder et al. 2009b). In T. jacksonii , re¬ 
production is clearly associated with higher precipitation 
during the long wet season in East Africa (Lin 1979, Lin 
& Nelson 1981). Although the amount of annual precip¬ 
itation is quite similar in the chameleons’ native and in¬ 
vaded ranges, precipitation seasonality is much more pro¬ 
nounced in its native range. Thus, the higher monthly pre¬ 
cipitation in its Hawaiian range may even facilitate circu- 
mannual breeding. 

Following the terminology of Swets (1988), we received 
good’ average AUC values for the SDM-based on current 
climate (AUC . . = 0.86; AUC + = 0.78). The minimum 
prediction value of the training records was 0.12, and the 
lowest 10 percentile was 0.27. Analysis of variable contri¬ 
butions in the model revealed that on average the precipi¬ 
tation in the coldest quarter’ had the highest explanative 
power (30.2 %), followed by the mean monthly temper¬ 
ature range’ (26.1 %), and the precipitation seasonality’ 
(12.8 %). Average contributions of all other variables were 
< 10 %. Our SDM suggests that, under current ‘bioclimat- 
ic’ conditions, Jackson’s chameleon may find suitable are¬ 
as nearly everywhere on the Hawaiian Islands, including 
many reserves, even when the availability of suitable mi¬ 
crohabitats in terms of land cover is considered (Fig. 3). 
The situation is not likely to change under future anthro¬ 


pogenic climate change with A2a and B2a scenarios, which 
forecast nearly identical and only slightly less suitable con¬ 
ditions for T. jacksonii (Fig. 4). 

Invasive Hawaiian Jackson’s chameleons are suggested to 
interact with native biota, at least as a predator to endemic 
fauna (Holland et al. 2009). While the impact of this spe¬ 
cies is poorly understood, we here provide evidence how 
severe the problem may be or may become as T. jacksonii 
finds suitable habitat almost all over the Hawaiian Islands. 
This emphasizes, in concert with Holland et al. (2009), 
the severity of a problem that should no longer be neglect¬ 
ed. Besides the impact on potential prey invertebrates, it 
should be considered hypothetically that T. jacksonii might 
compete with other local fauna for food resources. 

The potential invasion success and expected damage 
to the native fauna may be unparalleled by other invasive 
chameleons, although potential distribution and ecological 
impact studies remain limited (e.g. Dimaki et al. 2001). At 
least the spatial potential of Chamaeleo africanus in main¬ 
land Greece may be limited as, despite of the species’ oc¬ 
currence on the Peloponesos since historical times, there 
is no indication of past or present dispersal (W. Bohme, 
unpubl. data). Likewise, the Reunion population of Furci- 
fer pardalis is suggested to be restricted to a small, particu¬ 
larly dry and hot area on the western coast near St. Paul, 
although this species has recently been found elsewhere on 
the island, but its invasion success there remains unknown 
(E. Koch, pers. comm., March 2010). Invasion success was 
apparently limited in Furcifer oustaleti (Mocquard, 1894), 
which was translocated from Madagascar to near Nairobi, 
Kenya, but where this species has not been documented 
since 1974 (Spawls et al. 2002). Trioceros jacksonii seems 
to be exceptional among chameleons for its invasive poten¬ 
tial, as shown by its successful colonisation of the Hawai¬ 
ian Archipelago. 
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